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1. Introduction 

With the ongoing evolution of the experimental programs of the LHC and the Interna- 
tional Linear Collider, high precision predictions for multi-particle processes are urgently 
needed. In the last years we have seen a remarkable progress in the theoretical description 
of multi-particle processes at tree-order, thanks to very efficient recursive algorithms [1]. 
Nevertheless the current need of precision goes beyond tree order and therefore a similar 
description at the one loop level is more than desirable. 

The computation of the one-loop matrix elements seems to be notoriously difficult. 
The development of the main ingredients to accomplish this task, started almost 30 years 
ago with the pioneering works of 't Hooft and Veltman [2] and Passarino and Veltman [3] . 
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Still nowadays very few complete calculations with more than 3 particles [4] in the final 
state exist ^. 

The problem arises because of two reasons: the complexity of Feynman graph rep- 
resentation at the one-loop level and the way the reduction of n-point one-loop integrals 
in terms of scalar 4-, 3-, 2- and 1-point functions is performed. For the former, recur- 
sive equations seem to be a very promising tool. For the latter, it would be desirable to 
have a reduction of the full (sub-)amplitude instead of the individual tensor one-loop inte- 
grals. In fact, a method that will be based on the minimum possible analytical information 
about the one-loop amplitude will be more adequate, in principle, for an efficient numerical 
implementation^ . 

During the last years we have seen some very interesting developments. In the front 
of the reduction of tensor integrals a new method at the integrand level has been worked 
out by Pittau, and del Aguila [7]. This method will be the starting point of our work. 
On the other hand the idea of cut constructibility [8] in computing one-loop amplitudes 
(and not just integrals), has been proven very efficient in getting many very well known 
results. In the last two years a new development in computing one-loop amplitudes has 
been initiated by the work of Britto, Cachazo and Feng, [9] and combined with the unitarity 
method, resulted to a scries of remarkable results concerning QCD one-loop multi-particle 
amplitudes [10]. Moreover the introduction of qiiadruple cuts allowed a simplified algebraic 
approach, at least for the coefficient of the box function. 

In this paper we propose a reduction of arbitrary one-loop (sub-)amplitude at the 
integrand level by exploiting the set of kinematical equations for the integration momen- 
tum [11], that extend the already used quadruple, triple and double cuts. In contrast to 
the usual method of cut-contractibility, it is also possible to reconstruct the full rational 
terms of the amplitude. The method requires a minimal information about the form of the 
one-loop (sub-) amplitude and therefore it is well suited for a numerical implementation. It 
also gives rise to very interesting simplifications of well known results. 

In [7] it has been shown, by explicitly reconstructing denominators, that the integrand 
of any m-point one-loop amplitude can be rewritten as ^ 

A{q) = ^^^^ , Di = {q+pi)^-ml (1.1) 

where we use a bar to denote objects living in n = 4 -|- e dimensions and where the 
numerator N{q) can be cast in the form 

m— 1 m— 1 

^(q) = X] d{ioiii2i3) + d{q; iQiii^iz) JJ A 

io<ii<i2<i3 i^«0;«i,22,*3 
m— 1 m— 1 

+ X] [c(^on«2) + c(g;«on^2)] W Di 

io<ii<i2 ij^io,ii,i2 



^For a review of reduction methods, see [5]. 

^For other attempts towards a direct numerical implementation sec also [6]. 

^We assume po 7^ 0, for that choice allows a completely symmetric treatment of all denominators Di. 
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m—1 m—1 

+ J2 [K^on) + %;ioii)] n 

io<ii i^ioyii 
m—1 m—1 
+ [a{iQ)+a{q-io)] JJ A 

m—1 

+ P(g)nA. (1.2) 

i 

The quantities d(ioni2^3) are the coefficients of the 4-point loop functions with the four de- 
nominators Di^Di^Di^Di^. Analogously, the c(zo^ii2), 6(^0^1) and a(io) are the coefficients 
of all possible 3-point, 2-point and 1-point loop functions, respectively. 

The "spurious" terms d, c, a and P still depend on q. They are defined by the 
requirement that they should vanish upon integration over d^q, as we shall see later. 

Notice that no coefficient of scalar functions with more that four denominators appear. 
This is due to the fact that scalar functions with m > 4 can always be expressed as a linear 
combination of 4-point functions and, possibly, extra d terms [12]. 

All g's in the numerator N{q) of Eq. (1.2) are 4-dimensional. If n-dimensional t/'s are 
needed, they can be split into 4-dimensional and e-dimensional parts as explained in [7]. 
Eq. (1.2) is then applicable to the purely 4-dimensional terms, that are usually the most 
difficult to compute. 

Since the scalar 1-, 2-, 3-, 4-point functions are known, the only knowledge of the exis- 
tence of the decomposition of Eq. (1.2) allows one to reduce the problem of calculating A{q) 
to the algebraical problem of extracting all possible coefficients in Eq. (1.2) by computing 
N{q) a sufficient number of times, at different values of q, and then inverting the system. 

In carrying out this program two problems arise. First the explicit knowledge of the 
spurious terms is needed, secondly the size of the system should be kept manageable. To 
illustrate this second point, for m = 6 and without counting any spurious term, there are 
56 independent scalar loop functions and it is clearly advisable to avoid the inversion of a 
56 X 56 matrix. Our solution to this is basically singling out particular choices of q such 
that, systematically, 4, 3, 2 or 1 among all possible denominators Di vanishes. Then, as 
we shall see in Section 3, the system of equations becomes "triangular" : first one solves for 
all possible 4-point functions, then for the 3-point functions and so on. 

Notice that the described procedure can be performed at the amplitude level. One 
does not need to repeat the work for all Feynman diagrams, provided their sum is known. 
This circumstance is particularly appealing when our method is used together with some 
recursion relation to build up N{q). We postpone this problem to a future publication and, 
in this paper, we suppose to know N{q). 

A last comment is in order. In reconstructing the denominators, there is a mismatch 
between the 4-dimensional q in N{q) and the n-dimensional denominators D^. To compen- 
sate for this it suffices to replace m? mf — (f in all the coefficient of Eq. (1.2), where (f is 
the (n — 4)-dimensional part of (f' [7]. The coefficients of the various powers of q^, obtained 
through this mass shift, are the coefficients of the extra integrals introduced in [7], which 
give rise to nothing but the rational part of the amplitude [13]. 



-3- 



In this paper we show how to extract, with the help of Eq. (1.2), the coefficients of 
the loop functions, including the rational terms, from any amplitude. In Section 2 we list 
and compute the needed spurious terms, in Section 3 we show how to get, in a systematic 
fashion, the coefficients of all 4-, 3-, 2- and 1-point scalar integrals, while Section 4 deals 
with the problem of adding the missing rational terms. Finally, in Section 5, we present 
some practical applications and tests we made on our method. 

2. The spurious terms 

Before any attempt of extracting ci(ioiii2^3)) c(ioiii2)) K'^'oh) and a(io), one should ex- 
plicitly know the q dependence of the spurious terms d{q; 10111213) , c{q; ioiii2) , b{q;ioii), 
a{q;io) and P{q)- This can be achieved by decomposing any 4-dimensional q appearing in 
the numerator of Eq. (1.1) in terms of a convenient basis of massless 4-momenta [7], the 
coefficients of which either reconstruct denominators or vanish upon integration over d^g. 
The first terms gives rise to the d, c, b and a coefficients in Eq. (1.2). The second ones 
to all the additional spurious terms. The latter category is further classified, in Eq. (1.2), 
according to the number of the remaining denominators. In the following, we shall call the 
d, c, b, d, P terms 4,3,2, 1,0-point like spurious terms, respectively. As we shall see, the 
actual form of them generally depends on the maximum possible rank of the loop tensors 
in the amplitude. 

Since the decomposition makes use of the momenta appearing in the denominators 
we set, for simplicity, = 0, ii = 1, 12 = 2 and = 3 and derive explicit formulae 
for d{q;0123), c{q;012), b{q;01), d{q;0) and P{q)- By relabeling the indices, all the other 
spurious terms are easily derived. Before carrying out this program, we recall some basic 
results of Ref. [7], also with the aim to set up our notation. 
The explicit decomposition reads 



9' 




^[2[{q+Po)-ii]i^ + 2[{q + po)-i2]i1], 
[(9 + Po)■4]<+[(g + Po)•^4]^^ 



(2.1) 



where £1 and I2 are massless 4-vector satisfying the relations 



ki= £1 + ai£2 , k2= £2 + oi2£\ , 



(2.2) 



with 



ki=pi-po. 



(2.3) 



Furthermore, in spinorial notation. 



^^=<£i|7n^2], <=<^2|7n^i] with (4 • ^4) = -4(^1 • ^2) . 



(2.4) 
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The solution to Eq. (2.2) reads 

h = Pih - aik2) , h = /3{k2 - a2ki) , 

/3 = 1/(1 - aia2) , at = — , 

1 

7 = 2{h ■ h) = {ki-k2)±VA, A = {ki ■ k2f - klkl . (2.5) 

A last comment is in order. We make the assumption, always realized in practical 
calculations, to compute the amplitude A{q) in a gauge where the maximum rank of the 
appearing loop tensors is never higher than the number of denominators, as it happens, 
for example, in the rcnormalizablc gauge. This choice limits the number of the spurious 
terms. For instance, there is no P in such a case. However, the fact that all gauges are 
equivalent, leads us to the following conjecture 

One can always limit her/himself to the spurious terms appearing in the renoramaliz- 
able gauge, because, in physical amplitudes, the contributions coming from the tensors of 
higher rank should add up to zero. 

This conjecture, being rather strong, has to be checked in practical calculations. We are 
now ready to derive the q dependence of the spurious terms. 

2.1 The 4-point like spurious term 

By iteratively using Eq. (2.1), only one possible integrand that vanishes upon integration 
is left, namely 

d{q- 0123) = J(0123) T{q) , (2.6) 
where d(0123) is a constant (namely does not depend on q) and 

T{q) = Tr[i^ + Mhhhl5]. (2.7) 

To prove this statement, let us call N^^\q) the numerator of a term containing the four 
denominators DQD1D2D3. N^^\q) is necessarily a polynomial in q, whose maximum degree 
we will denote by jmax ■ Notice that jmax is also the maximum rank of the 4-point tensors 
that appear when performing a standard reduction procedure. Being interested in terms 
in which the four original denominators are not canceled out by the numerator function, 
we can systematically neglect all denominators that are reconstructed from N^^\q) ^. In 
particular, by expressing back £1 and £2 in terms of ki and k2, as shown in Eq. (2.5), one 
obtains [7] 

= |(po • + |(P0 • + {Di -Do + hi) r^ + {D2 - Dq + /is) , 

h^ = {m\ -p\)- {ml - vl) , = - ai£^ , r^^ = £% - . (2.8) 

*We will also ignore terms proportional to powers of q^, for they give rise to rational parts in the 
amplitude that can be treated separately (see Section 4). 
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Therefore 



1=0 



where 



so that 



2 



7 



(2.9) 



(2.10) 



(2.11) 



j=0 



We used the notation 0{Di) to indicate terms in which one of the denominators Di has 
been exphcitly reconstructed, and can therefore be neglected, as far as the construction 
of spurious terms is concerned. By replacing each q appearing in N^^^ (q) by the r.h.s. of 
Eq. (2.11), the only possible numerators of degree jmax, which preserve all four denomina- 
tors, turn out to be 

[{q + po) ■ isY-- , [{q + Po) ■ 14]'"- , [{q + Po) ' f^sY' [{q + Po) ■ hY' , (2.12) 

with j3 + j4 = jmax- The rank of such terms can be reduced with the help of the two 
following identities ^: 

[{q + Po) ■ 4] [(g + po) • h] = /3{q + PoTD^ - jiq + Po)' , 
[{q + Po) ■ h{4)][{q + po) ■ 4(4)] = { \rri<l + po)' - P{q + Po)"i5a] (ks ■ 4(4) ) 

- 2[7[(g + Po)-A:3]-/3A;^I?a][((Z + Po)-4(4)]} • (2-13) 
In fact, the insertion of Eqs. (2.9) and (2.11) in Eq. (2.13), together with the identities 

{q + pof = Do + ml-q\ 2 (g • fcg) = D3 - Do + hs , (2.14) 
with /13 given in Eq. (2.8), gives 

2^2 ^2 -pF Q 



[{q + Po) ■ h][iq + Po) ■ i^l 
[(q+Po) •4(4)][(« + Po) •%)] 



27 



1 



(^3 -4(3)) 

2 



(^3-4(4)) 



i=0 
32 Tp2 



2/3^ -PF-Q 
7 rng 



27 



[(g + po) -4(4)] 



+ YO{Di) + 0{f). 



(2.15) 



i=0 



■'a demonstration can be found in [7]. Notice that the factor 2 in front of the last term of the second 
equation is missing in that paper. 
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The two equations (2.15), when introduced in the numerators of Eq. (2.12), reduce their 
rank from jmax to jmax~^j up to contributions containing less denominators or proportional 
to q^. By applying the described procedure jmax — 1 times N^^\q) takes the form 

3 

iV(3) {q) =r,o + m [{q + po) ■ 4] + m [{q + Po) " 4] + ^( A) + 0{f) , (2.16) 

i=0 

where the coefficients riQ, r/3 and 774 do not depend on q. The denominators still hidden in 
[{q +Po) ■ ^3,(4)] can be further extracted by using the identity [7] 

[{q + Po) ■ 4,4] = (^^^ [pktD^ -l[{q+ Po) " ^s] ± ^ } 

pk3-F--f 



1 



Po-h + Y 



(^3 • 4,3 

3 

+ Y,OiDi) + 0{e), (2.17) 

i=0 

with T{q) given in Eq. (2.7). Therefore the final expression for N^^\q) reads 

3 

N^^\q) = d(0123) + d(0123) T{q) + ^ 0(A) + 0{f) . (2.18) 

i=0 

The statement that d{q; 0123) must have the form given in Eq. (2.6) is then equivalent to 
the 

Theorem : 

d"q -^^i- = . (2.19) 

The proof trivially follows by making the shift q ^ q — po in the integration variable and 
by noticing that T{q) oc e{q,li,i2,ks). In fact, the resulting rank one 4-point function can 
only be proportional to fci, ^2 and ^3 and each term vanishes when contracted with the e 
tensor. 

2.2 The 3-point like spurious terms 

To derive c{q; 012) we make use of the 
Theorems : 

Once again, they are proven by shifting q and performing a tensor decomposition. 

Such integrals come from the only possible terms in the decomposition of 3-point 
tensor integrands that do not reconstruct denominators [7]. The proof closely follows the 
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reasoning of the 4-point case. Calling N^'^\q) the numerator of a term containing the 
three denominators D0D1D2, and being jmax its maximum rank, one applies Eqs. (2.8)- 
(2.11) to cast N(~'^\q) in a form where all rank jmax 

terms with three denominators are 
proportional to the three numerators given in Eq. (2.12). The first two terms are not 
further reducible, while the iterative use of the first of Eqs. (2.15) reduces the third one 
to a sum of contributions proportional to [{q + Po) •^3]-' and [{q + Po) ■ separately, with 
j < jmax- Then 

Jmax 

iV(2)(g) = c(012) + {ci,(012)[(g +po) • + C2j {012)[{q + po) ■ £4]'} 
2 

+ YO{Di) + 0{f), (2.21) 

i=0 

where c(012), cij(012) and C2j(012) are constants. Therefore c(g; 012) must have the form 

Jmax 

c{q; 012) = ^ {cij (012) [{q + po) ■ + C2j (012) [{q + po) ■ 4]^' } ■ (2.22) 

In the renormalizable gauge, jmax = 3. In order to illustrate this fact, a simple argument 
can be used. Let us consider the reduction of a m-point function of rank m. In the 
decomposition of Eq. (1.2), the 3-point like spurious terms involve (m — 3) reconstructed 
denominators -Dj. Each one of them can be obtained from a power of q in the numerator 
N{q) by means for example of Eq. (2.11). This leaves at most m — [ra — 3) = 3 powers 
of q available for the construction of the c{q) term. The same reasoning can be applied to 
determine jmax = 2 and jmax = 1 respectively for 2-point and 1-point like spurious terms. 

2.3 The 2-point like spurious terms 

To derive h{q\ 01) it is convenient to rewrite g in a basis expressed in terms of an auxiliary 
arbitrary 4- vector v not parallel to ki ^: 

= -P^ + yikt + Vun" + yjt^ + ■ (2-23) 

In a way similar to Eqs. (2.2) and (2.4), ki and v are decomposed in terms of two massless 
4- vectors I^^q, 

ki= £5 + "54 , v = £e + adb , (2.24) 

and £7,8 defined as follows 

^^ = <4|7l4], £^=<£6\r\£5]. (2.25) 

Then n is taken to be 

n = 4 - "54 , (2.26) 

® We normalize w such that (ir-is) = -4(4-4) = -kj. Then as = 2, £5 = (fci+n)/2 and 4 = (fei-n)/4. 
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so that it satisfies the two following properties 

n-ki=0, and = -kf . (2.27) 

Finally, one computes 

, [{q + Po) ■ ki] ,^1, [{q + Po)-n] ^^ , [(g+Po) -4] ,/. , [{q+Po)-h] ,^ 

(2.28) 

By using this basis, the spurious terms can be determined with the help of the following 
Theorems : 

^n J(g+P0)-^7p'[(g + P0)-n]' _ 

^ DoD^ 
^- [{q + Po)-h]niQ + Po)-nY ^ 
DqDi 



I 



I 



I 



,n- [iQ + Po)-n?'-' 



[{q+po) • n]^^ - rj{[{q+po) • fci]^ - (g +Po)^fcf p 11 

Vj = l,2,3,--- and i = 0,l,2--- (2.29) 



Again shifting q and decomposing allows an easy proof. As for the fourth line of Eq. (2.29), 
both terms in the numerator, after tensor decomposition, are proportional to k^-^ . The 
factors rj are chosen such that their STim is zero. 

Now we are ready to prove that only the terms given in Eq. (2.29) contribute to b{q; 01). 
First one rewrites 

= G''-T2 {[{<!+ Po) • njn'^ + [{q + Po) ■ ir^ + Hi + Po) ■ 4]^^} + J] 0( A) , 



kf 



(po • ^i) + Y 



i=0 

(2.30) 



with hi given in Eq. (2.8). Then, calling N^^\q) the numerator of a term containing the 
two denominators DqDi, and being jmax its maximum rank, one replaces each q appearing 
in N^^\q) by the r.h.s. of Eq. (2.30). Therefore, the only possible generated numerators 
of degree jmax are 

[{q + Po)-i7y, [(g+Po)-4P, 

[{q +po)- iiV' [(q +Po)- nY- , [{q+Po)' [(g + Po) • n^" , (2-31) 
[{q + Po)- ny, [{q + Po) ■ ^iX' [{q + Po)- kY' , 

with j = jmax, j7,8+j n — jmax and + ig — jmax- Duc to the first two lines of Eq. (2.29), 
the first four terms of Eq. (2.31) directly give rise to 2-point like spurious terms. Conversely, 
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the last one can be further reduced by means of the following identity 



[{q +Po)- £7] [{q + Po) • 4] = 4[(g + po) ■ h] [{q + Po) ■ 4] - 2(4 • 4) {q + Pof 

= ^ I (iPo • fci) + y) ' - kWo - [{q + Po) • 1 

1 



(2.32) 



i=0 



where we have used the results of footnote 6 and reconstructed denominators as in Eq. (2.14). 
Then, the last term of Eq. (2.31) can be put in the form 



[{q + Po) ■ irr [{q +Po)-isr = Y1 c'( A) + o{f) 

) 

[{q+Po)-i8Y'-''Y^Si[{q+po) 



+ < 



«7 



n 



l2i 



i=Q 

is 



[iq+Po)-i7Y'-''^Si[{q + 



Po)-nf, 



(2.33) 

if 17 < is, 
if ^8 < ir, 



i=0 



where the 5i are constants. All pieces generated by the previous equation are given by the 
first five terms of Eq. (2.31), but now with j < jmax and jy.s + in ^ jmax- Therefore, apart 
from [{q + po) ■ n]^ , they are again taken into account by the first two lines of Eq. (2.29); 
in other words, they contribute to b{q; 01). If j is odd, due to the third line of Eq. (2.29), 
[(q +Po) ■ also gives rise to spurious 2-point terms. On the contrary, in order to get the 
contribution to b{q; 01) in the case when j is even, one has first to subtract powers of 



[{q+Po)-h\'-{q+pofkl 



(2.34) 



as performed in the last line of Eq. (2.29). Incidentally, when adding this piece back, a 
contribution to 6(01) is generated, since 



[{q+PQ)-h]''-{q + pofkf 



{po ■ h) + 



— kfrriQ 



i=0 



(2.35) 



In summary, no other 2-point like spurious terms are possible besides those listed in 
Eq. (2.29). 

In the renormalizable gauge there are at most two powers of q in the numerator. By 
applying the above reasonings with jmax = 2 and jmax = 1 gives rise to eight possible 
spurious b terms: 

6(g;01) = 6ii(01)[((?+po)-4] + &2i(01)[(g+Po)-4] 
+ 612(01) [((? +po) • 4]' + b22i01)[{q + po) ■ is? 
+ 6o(01)[(g + Po) ■ n] + 6oo(01) Kiq; 01) 
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+ boiiOmq + po) ■ i7][{q + po) ■ n] 

+ 6o2(01)[(g +Po) • + Po) • n] , with 



i^(.;01) = { [(q + po) • - KilMl^ll^ziil^ ^ . (2.36) 



2.4 The 1-point like spurious terms 

First we decompose 

q^^ = -p^ + yk^' + y^ni" + y7ii} + ysi^ , (2.37) 

where k is an arbitrary 4- vector and n, ij and is are built up from k and v as in the 2-point 

case. Then we make use of the following 

Theorems: 



/ 



Do 

Do 

Tn = (5miM2-M2„5'^^''' • • • g^^'-'^'-r^ , Vn = 1, 2, 3, • • • , (2.38) 

for any 4-vector Vi and where 5/ii/i2 - M2n the symmetrized product of n metric tensors. 
The proof is just a direct consequence of the fact that 



/ 



(^Q^—T^ = and d^q^—^ 1^ oc 5'^i'^2 -'.2„ (2.39) 

(r - ml) J [q^ - ml) 



In the renormalizable gauge at most rank one 1-point functions appear. Therefore 

a(g;0) = diiO)[iq + Po) ■ k] + d2i0)[{q + Po) ■ n] 

+ a3(0)[(g+Po)-^7]+a4(0)[(g+po)-4]- (2.40) 

2.5 The 0-point like spurious term 

As already pointed out, P{q) = in the renormalizable gauge. We can prove this statement 
simply by counting the powers of q in Eq. (1.2). The last term on the r.h.s. contains 2m 
powers of q in the m reconstructed denominators Di. Since N{q) on the l.h.s. is at most 
of rank m and the other terms on the r.h.s. contain at most 2m — 1 powers of q, in order 
to satisfy Eq. (1.2) we should have P{q) = 0. 

In more general gauges, P{q) gives a contribution, polynomial in q, to the integrand 
of the amplitude. Terms like that vanish, upon integration, in dimensional regularization. 
This is why we classified P{q) among the spurious term. Off course P{q) is never needed, in 
any gauge, «/all the other coefficients in Eq. (1.2) can be determined without making use of 
its actual form. This is the case, as we shall see in the next Section. However, if necessary, 
after all the other coefficients are known, P{q) is easily computed as the difference between 
N(q) and the first four lines of the r.h.s. of Eq. (1.2), divided by all the denominators. 
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3. Extracting the coefficients of the scalar loop functions 

Being interested here in the coefficients of the scalar loop functions, we can set everywhere 

g' = 0, (3.1) 

so that, in particular 

Di^Di = {q + pif-mj. (3.2) 

The "error" induced by the above replacement is at the level of the rational part of the 
amplitude, as we shall see in Section 4, where we will also use the fact that the terms 
are always connected to the masses in the denominators to reconstruct the information we 
are missing with the replacement of Eq. (3.2). 

Our approach to the problem is suggested by the structure of Eq. (1.2) itself. Choosing 
particular values of q such that 4, 3, 2 or 1 denominators vanish allows one to reduce the 
problem to the solution of simpler sub-systems of equations. To illustrate how this works 
we concentrate again on the particular choice io = 0, ii = 1, ^2 = 2 and is = 3 and derive 
explicit formulae for the coefficients d(0123), c(012), 6(01) and a(0) as well as for the 
coefHcients of the corresponding spurious terms ^. As we shall see, the latter information 
is also needed when iterating the algorithm: in order to solve for the c(012) coefficient, one 
has to know the coefficients of all terms with 4 denominators, including the spurious ones. 
In order to solve for 6(01), one needs, in addition, all the terms with 3 denominators and 
so on. 

3.1 The coefficient of the 4-point functions 

We look for a q such that 

Do = Di = D2 = Ds = 0. (3.3) 

By writing the 4-vector q as 

4 

q^ = -Pl; + J2^^i';, (3.4) 
1=1 

with £i given in Eqs. (2.2)-(2.5), one obtains a system of equations for the xf. 

= "f{xiX2 — 4x3x4) - do 
= do ~ di + 7(.Tiai + X2) 
= do - d2 + 7(2^202 + xi) 

= do - da + 2 [xiih ■ h) + 0:2(^3 • ^2) + x^ih ■ 4) + 2:4(^3 • 4)] , (3.5) 
where ki is given in Eq. (2.3) and 

d, = mf - A;f . (3.6) 
'^For any other choice of indices the procedure is the same. 
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There are two possible solutions 

{qo^)i^ = -p^ + x'ie^ + x^^i^ + X3^e^ + X4^i'l , (3.7) 



with 



Xi = —[d2 - a2di - do{l - a2)] ■ 
7 

X2 = —[di — aid2 — do{l — ai)] - 

7 

Axf^ + Bxf -C = 0, 

±_ c 
x^ 



and 



X3± 



, _ jks ■ 4) ^3 - do - 2x°(fc3 • li) - 2x°(fc3 • i2) 



(3.8) 



Notice that we need two solutions to be able to determine both d(0123) and (i(0123), and 
that the existence of more than one solution is a consequence of the quadratic nature of the 
system in Eq. (3.5). By putting both in Eq. (1.2), and recalling the form of d{q; 0123) 
given in Eq. (2.6), one finds 

iV(go±) = [d(0123)+J(0123)r(go±)] [] Di{q^) . (3.10) 

17^0,1,2,3 



Then, by defining 



it is possible to extract d and d 

,(0123) = »^41^^^^£1^, 

d(0123) = ||4l^|M. (3.12) 

Notice that, in terms of 0:^4, one rewrites 

T{q^) = 27 [xf{k3 ■ 4) - xi{k3 ■ £4)] , 

T{q+) = -T{q^). (3.13) 

Then 

d(0123) = ^[R{q+)+R{q^)], 

m2S)='-^^^^l^. (3.14) 



-13- 



The two above equations do not depend on the rank of the tensors in the amphtudc. When 
N(ci) = 1, they allow a trivial decomposition of any m-point scalar loop function with 
m > 4 to boxes, as we shall see in Section 5. 

3.2 The coefficient of the 3-point functions 

At this stage all d and d coefficients are known. When q is such that 

Do = Di=D2 = Q and A 7^ 7^ 0, 1, 2 (3.15) 

Eq. (1.2) reads 

N{q) - [^(0i2i3) + 0i2i3)] n ^^('?) 

2<i3 1^^0,1,2,13 

= R'{q) n A(g) = [c(012)+c(g;012)] [J A(g) , (3.16) 

17^0,1,2 i^0,l,2 

and one can extract c(012), together with all the six Cij(012) coefficients of Eq. (2.22), by 
computing R'{q) at seven different q's that fulfill Eq. (3.15). For a q written as in Eq. (3.4) 
that happens when 

xi = x\ 

X2 = X2 

X3X4 = C , (3-17) 

where X12 s^nd C are given in Eqs. (3.8) and (3.9). There is now an infinite number of 
solutions, which we parametrize by imposing the extra condition 

{q+Po)-h = ±VCe'''/''{i3-e4), k = 1,2,3,- ■■. (3.18) 

Then 

X4 = ±VCe'^/'' = x^f^, X3 = -^ = ±VCe-'^/^ = x^^, (3.19) 

and the complete solution reads 

(Q^r = -P^o + ^Ti + ^^^^ + 4/1 ■ (3.20) 
Finally, with the help of the relation 

[{qt + Po) ■ ^4] [{q^ + Po) ■ 4] = C(4 • , (3.21) 

one writes 

R'ili)= E c,[±e-/Y, (3.22) 

with 

In Appendix A, we explicitly determine all Cj's of Eq. (3.22) by choosing 7 different solutions 
of the form given in Eq. (3.20). 
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3.3 The coefficient of the 2-point functions 

At this stage all d, d, c and c coefficients are known. When q is such that 

Do = Di = and A 7^ Vi 7^ 0, 1 (3.24) 

Eq. (1.2) reads 

N{q)- Yl [d{01i2i3)+d{q; 011213)] J] 

I<i2<j3 1^0,1,12,13 

- 5^[c(01i2)+c(^;01i2)] n 

l<i2 1^0,1,12 

= R"iq) n Di{q) = [6(01) + b{q; 01)] J] A(g) , (3.25) 

and one can extract 6(01) together with all the eight 6(01) coefficients of Eq. (2.36), by 
computing R"{q) at nine different g's that fulfill Eq. (3.24). That happens when, for a q 
written as in Eq. (2.23), 

di — do Q 

vl^ty^f-"!^^^, (3.26) 

where we used our normalization condition {£7 ■ £§) = —kf. We impose now, as an extra 
requirement, that K{q; 01) defined in Eq. (2.36) vanishes . This implies 

yi = yi, 



2/7 = - • (3.27) 

2/8 

We fix the remaining freedom by imposing 

{q + Po) ■ h = ±VFe'''/''{£7 ■ 4) , (3.28) 

that implies 

2/7 = ^ = iVFe-^'^/*^ = y^fc . (3.29) 

Then, this class of solutions can be parametrized as follows 

& = -P^ + yikt + l^n^ + y^k ^7 + ytk ^8 > with l = ±l. (3.30) 
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Since 



W + Po)-4 = ^— — 3.31 



we can rewrite (see Eq. (2.36)) 



R"{ifk) = iY.^3 [^^'^^^y + E ^3 [^^'^^^y , (3-32) 

j=-l j=-2 



where 



6(01) = 60 bijm= , 62,(01) ^" 



^oi(Ol) = ^tt/W bo2{01) = -J^^ bom ^° 



(3.33) 



To determine the last coefficient 600 (01) we need to introduce a different solution, for which 
K(q; 01) does not vanish. We call qq such a solution and choose it satisfying the condition 
2/7 = 2/8 = 0. Then, by insertion in Eqs. (3.26) and (2.36), one finds 

q^ = -p^ + y^ik'^ + V3Fn'', (3.34) 

and 

R"{qo) = 6(01) + ^o(Ol) [(go + Po) ■ n] + 6oo(01) K{qo; 01) 

= 6(01) - 60(01) VsFkj + 2600 (01) kjF . (3.35) 

The complete solution for all the coefficients is given in Appendix B. 

3.4 The coefficient of the 1-point functions 

In massless theories all 1-point functions, namely all tadpoles, vanish, also implying that, 
in such cases, one does not need to know all the 6 coefficients: only the coefficients of the 
scalar 2-point functions are needed, each of them can be determined in terms of just four 
solutions of Eq. (3.24) (see Appendix B). However, in general, also the coefficients of the 
tadpoles are required. Therefore we show how to extract them. 

At this stage we assume to know all the d, d, c, c, 6 and 6 coefficients and, when q is 
such that 

Do = and A 7^ Vz 7^ , (3.36) 

Eq. (1.2) reads 

N{q) - Yl [d{0iii2i3)+d{q;0iii2h)] JJ 

0<ji<fi2<j3 iy^Q,ii,i2,i3 

- [c(0iii2) + %;0zii2)] n 

0<ii<i2 1^0,11,12 
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0<-ii ij^O,ii 

^ R!"{q) n = [a(0) + ~a{q- 0)] J] A(«) • 



(3.37) 



Then, we parametrize q as in Eq. (2.37) and choose two solutions of Eq. (3.36) such that 

yn = y7 = ys = O: 



{q^Y = -i^o^\l^k^^. 



Therefore, by recalling Eq. (2.40), we can write 



El"{q+) = a{Q) + am\^k' 



fc2 

B!"{q^)=a{Q)-~am\i%k\ 



so that 



a(0) = 



B!"{qt) + B!"{q^ 



(3.38) 



(3.39) 



(3.40) 



Notice that a(g;0) is never needed, because it would be necessary only to extract -P(q'), 
that, as already observed, is irrelevant. 

4. Reconstructing the rational part of the amplitude 

Until now we have assumed = 0. As already discussed, this is enough to reconstruct 
the coefficients of the 4-3-2- 1-point loop functions, but rational parts are missing. In the 
renormalizable gauge the only possible contributions to those rational terms come from the 
following extra scalar integrals introduced in [7] ^ 



/ 



dTq 



DiD.Dk 



I 



(fq 



DiDj 



ITT 



+ 0{e) 



+ rrij — 



2 (Pi-Pj] 



+ 0{e) . 



(4.1) 



We checked that they reproduce the rational terms listed in [13]. Therefore, in our language, 
the coefficients of the integrals in Eq. (4.1) are just the coefficients of the maximum powers 
of q^ contained in the d{ijkl), c{ijk) and b{ij) once q^ is reintroduced through the mass 
shift 



rrii 



~2 



(4.2) 



*The powers of are dictated by the maximum rank of the loop tensors in A{q). 
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With the above replacement the coefficients get a dependence on q and one can expand: 

d{ijkl;f) = d{ijkl) + fS^\ijkl) + qU^^\ijkl) , 
c{ijk;f) = c{ijk) +fc^'^\ijk) , 
biij;f) = biij)+fb^''Hij)- (4.3) 

d^^\ijkl), c^'^^ijk) and b^'^\ij) are then the coefficients of the first, second and third 
integral of Eq. (4.1), respectively. They can be either computed numerically 



d^^\ijkl) = lim 



d{ijkl; I 



c^^^ijk) = hm ^ 



b^%j)= hm (4.4) 
g^— >oo q 

or as solutions of systems obtained by evaluating Eq. (4.3) at different q^. For instance: 

n d{ijkl; 1) + d{ijkl; -1) - 2d{ijkl) 
d^ >{ijkl) = ^ , 

c(2) {ijk) = c{ijk; 1) - c{ijk) , 
b^''\ij) = b{ij;l)-biij). (4.5) 

A small example of the described approach is given in the next Section. 
5. Applications and tests 

We tested the whole method on the reduction of a rank four 4-point tensor integral 



/ 



an^JXtSL (5.1) 



to scalar functions. We have been able to correctly extract the coefficients of the scalar 
integrals. In addition, when reducing the tensor in Eq. (5.1) with the techniques of [7], we 
have been able to also test the coefficients of the spurious terms. 

As for the rational terms, we explicitly show here, as an illustrative example, the 
extraction of the coefficient of 



D0D1D2 



(5.2) 



from 



where, for simplicity, we have put po = 0- Then, in this case 

N{q) = q^q" = R'{q) . (5.4) 
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We first observe that, in the hmit — ^ oo, one can use the fohowing asymptotic form of 
the solutions given in Eq. (3.20) 



with Coo = lim C 

>oo 



4j 



Therefore, one obtains, for the coefficient cq in Eq. (A.4) 

Co _ ^3^4 + ^4^3 



Then 



.(2) 



(012) 



47 



_ ^3^1.4 



47 



(5.5) 



(5.6) 



(5.7) 



in agreement with the result obtained in [7]. 

Another rather straightforward application of the method is the reduction of the scalar 
n-point functions, with n > 5, in terms of box functions. It will allow the reader to follow 
the reasoning of the reduction in a simple case. It should be mentioned that the content of 
this derivation, amazingly enough, goes back to the year 1965, in the work of Melrose [14] 
and Kallen and Toll [11]. 

In the conventional approach one can easily prove [14, 15] the following relation 



(5.8) 



jN 


-7^-1(0) 


-7^-1(1) .. 


. -I^-\N-1 


1 


Yoo 


Yoi .. 


^OAf-l 


1 


Yio 


Yii .. 


YiN-l 


1 


Yn -10 


Yn-11 .. 


Yn-in-i 



where 7^ is the iV-point scalar function 



7 



N 



1 



'Do---Dn-i ' 

is the N — 1-point function with the i-th propagator missing and 
Yij = m^ + m]-{pi-pjf, i = 0,...,N j = 0,...,N. 



(5.9) 



(5.10) 



By repeated use of Eq. (5.8) we may express the A/"-point function in terms of 4-point 
functions with coefficients expressible in terms of the determinants of Y matrices. For 
instance 

'■^^^^7^(0, (5.11) 

t=0 ^ ' 



where Y^'^^ is the 5x51" matrix, and deti represents the determinant of matrix Y where 
all elements of the i-th column have been replaced by 1. 
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Similarly we obtain 

5 

det{Y^^)) 



i=0 



Let us now see how these formula get simplified using the method described so far. 
For the 5-point function we get 

4 

l'> = Y,dil\i) (5.13) 

i=0 

with 

whereas the 6-point function reads 

i<j i,j=0 

where is obtained from I® by dropping the propagators Di and Dj and where 

= 2 (a(4,))^.(4)) + mM^m^) • ^'-''^ 

In the above equations q^.-^ arc the two solutions given in Eq. (3.7) when all the propagators, 
except Di, are zero. Analogously q^^^ are the solutions when Di and Dj are the only non 
vanishing propagators. 

It is quite straightforward to prove that 

l(^_ l_\^_deU{Y(^ 

2 l^A(g(^)) D,iq-^)) det{Y(^)) ' ^ " ^ 

This is because the Gram determinant of q,Pi,P2,P3,P4,P5 should be zero, which results 
to 



2D0 + Y00 Di-Do + Y^o-Yoo D^-Do + Y^o-Yo^ 

Di — Dq + YiQ — Yqq Yii — Yio — Yqi + Yqo Y12 — lio — ^02 + ^0 
D2 — Do + Y2Q — Yqq Y21 — Y20 — Yqi + Yqo Y22 — Y20 — Y02 + Yq 



'00 • • • -Ds — Dq + Y20 — Yqq 
00 • • • ^15 ~ ^10 ~ ^05 + ^00 

...Y25- Y20 - Yo5 + Yoo 



D5-D0+Y20-Y00 Ysi-Y^o-Yoi+Yoo Y,2-Yso-Yo2 + Yoo . . . Y,^ -Y,o -Y05 + Y, 



00 



(5.18) 



when on takes into account that 



2q^ = 2Do + 2ml 
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2q-pj = Dj-Do + Yoj-Yoo, 

and 

2pi ■ pj = Yij - Yio - Yoj + Yqo . 
Taking Eq. (5.18) at the point q = q'(j) we end up with a second order equation for Dj given 

by 

aDf + 6 A + c = , 

with b = —2deti{Y^^^) and c = det{Y^^^). An analytical proof for arbitrary N, can be 
found in [14]. 

6. Conclusions 

We have shown how computing the integrand of any one-loop amplitude at special values 
of the integration momentum allows the one-shot reconstruction of all the coefficients of 
the scalar loop functions and of the rational terms. Then, by simply multiplying those 
coefficients by the known scalar integrals, the computation of the amplitude becomes trivial. 
Our method should be particularly useful in the case when recursive techniques are used to 
numerically compute the integrand. We plan to investigate this subject in the near future. 
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Appendices 

A. The system for the coefficients of the 3-point functions 

We choose the following seven solutions (see Eq. (3.20)) 

qf ^ it ^ it and , 

and define the combinations 

Then, from Eq. (3.22), one obtains, for the even powers of the system 

^^(^l) = +C-2 +C0 +C2 

< T+{q2) = -C_2 +Co -C2 , 

^ T+{q^) = +c_2e-2^'^/3 +C0 +C2e2^^/3 
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(A.l) 
(A.2) 

(A.3) 



whose solution is 
Co 
c±2 



T+(gi)+T+(g2) 
2 

r+(gi)-r+(g2) 



g±2»V3(r+(g3)-Co) 



\ _ g=F2i7r/3 



For the odd powers of j one gets instead 



T-iqs) 
T-{q2) 
T-{qi) 



-c_3 +c_ie-^''/3 +cie*''/3 

+iC_3 -zc_i +ici -iC3 

-C-3 -C_i -Ci -C3 

-zc_3 +c_ie-'''/6 +cie*''/6 +ic3 



where 



(A.4) 



(A.5) 



(A.6) 



is known because the coefficients co,±2 have aheady been determined. The solution reads 
[T-(g3) - r-(gi)](l + e±-/3) ^ i[rO(^g) + r-(g2)](l + e^-/3) 



c±i 



C±3 = -Cipl 



r-(gi)Tir-(g2) 



B. The system for the coefficients of the 2-point functions 

We choose the following nine solutions (see Eq. (3.30)) 

and build up the combinations 

Therefore the /? coefficients satisfy the system 



5- 


{qt) = 


-P-i 


+Po 


-Pi 




(Qt) = 


-if3-i 




+iPi 


s- 


(Qt) = 




+Po 


Pie'"" 



whose solution reads 

(1 + e^-/3)[S-(g+) - S-{qt)] - (1 T i)[S-{qt) - S-{qt 



P±i = 

Po = S-{qt) + {P-i+Pi) 



±i{3 - V3) 



(A.7) 



(B.l) 



(B.2) 



(B.3) 



(B.4) 
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Next, by defining 



one finds 



whose solution is 



T-{qi) = -bi 
T~{q2) = -ih-i +ibi 



(B.6) 



b±, = -l[T-{q^)±iT-{q2)] . (B.7) 



2 

Now that Po, 6_i and bi are known we define 

1 

T\q^) ^ R"iqX,^) - ^ (5j [e^^/^Y - 6_i e— - 61 e^/^ , (B.8) 



then 



= +6-2 +60 +^^2 

r+(g2) = -b-2 +bo -b2 ■ (B.9) 

This system is analogous to the system in Eq. (A.3) with the replacements Cj — and 
^^(93) T^{Q3)- Therefore its solution can be directly read from Eq. (A.4): 

r+(gi) + r+(g2) 



bo = 
b±2 



\ _ g=F2i7r/3 



(B.IO) 



Finally, from Eq. (3.35) one obtains the last coefficient 



Soo(oi)^ ^"<^°'-^;;-^^° . (B.n) 
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